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We analyse the production of a Higgs boson in association with a Z boson at hadron colliders in 
the Standard Model and some simple extensions. We show how multi-jet merging algorithms at 
leading and next-to-leading order for the loop-induced gluon fusion and the Drell-Yan like quark- 
induced processes, respectively, improve the descriptions for various differential distributions, in 
particular those that involve the production of additional jets. The phenomenological studies focus 
on two relevant channels of Higgs boson decays, namely H —>■ invisible and H —>■ bb. We find 
sizable and phenomenologically relevant corrections to the transverse momentum and invariant 
mass distributions for the Higgs boson candidate. Thanks to the large destructive interference for 
the top Yukawa terms, this process is very sensitive to the magnitude and sign of a possible non¬ 
standard top-Higgs coupling. We analyse the impact of this anomalous interaction on distributions 
and estimate constraints from LHC Run II. 


I. INTRODUCTION 

The preeminent achievement of Run I of the LHC was the discovery of a scalar particle resonance [1] , which so 
far proved largely consistent with the Standard Model (SM) Higgs boson [5113]. This discovery not only marks the 
end of an era of searches for this elusive resonance, but it also heralds the beginning of a new era of exploration 
of the electroweak symmetry-breaking mechanism. The increased collision energy and luminosity of LHC during 
Run H allows, in particular, precise measurements of the interactions of this new resonance with other known 
particles. At the same time, other new resonances interacting with the rest of the SM through the Higgs boson, 
and new structures in the interactions of known particles will become a primary ground of renewed rigorous 
searches for physics beyond the Standard Model. 

In this scenario, the associated production of a Higgs boson with a Z vector boson, pp ZH, also known 
as Higgs-Strahlung, is one of the most prominent paths towards an accurate understanding of the Higgs boson 
couplings. Remarkably, this production mode supplemented by jet substructure techniques can help to access the 
largest yet most challenging Higgs decay channel —>■ [4], whereas the leading gluon fusion and vector boson 
fusion channels fail in this task due to overwhelmingly large QCD backgrounds. ATLAS and CMS already have 
reported first hints for this process 13 HI; while the former collaboration provided an upper limit on the event 
rate of 1.4 times the SM expectation, the latter observed an excess of events above the SM background with 
2.1(7. Run H will thus clarify the situation concerning this process, fully establishing its existence and scrutinising 
its dynamics. Ultimately, the ZH channel will shed light on the highly relevant branching ratio of Higgs bosons 
decaying into invisible final states, an important portal for interactions between the Standard Model and the Dark 
Matter sector. This channel provides one of the strongest constraints, where the current upper bounds at 95% 
CL reported by ATLAS and CMS are BTZ{H —>• inv) < 0.75 and 0.58 [7], respectively. 



Figure 1: Feynman diagrams for ZH production at leading order: (a) Drell-Yan-like; (b,c) gluon fusion. 

In the SM, ZH production is dominated by the Drell-Yan-like mode, see Fig. [^a). At leading order (LO), it 
contributes to the total cross-section at Another relevant production mode of ZH final states is gluon 

fusion, a loop-induced process mediated by quark loops, depicted in Figs. [^b,c) and contributing at LO with 
0(a^a^^r). These contributions have been discussed for example in |8]. 

At the level of the Feynman diagrams shown, these two sub-processes do not interfere, but it is important to 
stress that the latter, gg —i ZH, is part of the next-to-next-to-leading order corrections to the total ZH production 
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cross section. Here, we will treat the two process classes as separate categories, since this allows a more careful 
study of the respective QCD emission patterns. There are four major factors that guarantee that the gluon fusion 
process is larger than the anticipated naive « 1%: i) it has a larger initial state colour factor; ii) the process 
is driven by the large gluon parton distribution function (PDF); iii) the top Yukawa coupling appearing in 

the box diagram in the place of one of the asw factors, is of order unity ~ C>(1); and iv) the top-quark 
loop presents a threshold enhancement at mzH ~ Sm*, which gives rise to relevant rates at the boosted regime 
Pth ^ mt- 

On the phenomenological side, and in particular in the framework of Higgs boson coupling fits, the loop-induced 
contribution provides an additional probe to the size and the sign of the top-quark Yukawa coupling. In Figs.[^b,c) 
the respective Higgs boson vertices lead to linear terms in Kt and Ky, where Kt is a New Physics deviation to the 
top-quark Yukawa coupling and Ky represents a potential rescaling of the HVV interaction vertices 

yt = gHvv = Kvgfivy > (1) 

with V = Z, W. On the other hand, the WH and qq —)■ ZR processes probe a single coupling strength, kw and kz, 
respectively. At the LHC there are other known processes able to probe both their size and sign, e.g., Higgs boson 
production in association with a single top quark pp —>■ tHj and the off-shell iL —>■ 4 lepton production [SlinHn]. 
However, their experimental observation is challenging due to small rates and huge backgrounds. 

In many analyses with complex final states, such as the ones emerging from the ZH processes discussed here, 
it has become customary to consider signals and backgrounds in bins of jet multiplicities; one of the most obvious 
examples being the process H —>■ WW, where the dominant tt WWbb background can be fought with jet 
vetoes. Similarly, albeit less importantly, the same logic can also be applied to the Rl~bb final states typical for 
ZH production. There, jet vetoes can play a role similar to considering boosted topologies, which also suppress 
the ti and similar backgrounds. This motivates a more detailed study of jet emission patterns in this process, 
where the tool of choice is the multi-jet merging technology that has already been used in a large number of Run I 
analyses, based on leading-order matrix element calculations. 

In this work we study improvements arising from multi-jet merging techniques applied to the simulation of Higgs- 
Strahlung, and the impact of the improvement which these techniques have experienced through the inclusion of 
next-to-leading order accurate matrix elements. The simulation comprises the following contributions: 

• The Drell-Yan-like pp HZ{ll) -I- 0,1 jets at next-to-leading order (NLO) accuracy in QCD merged into a 
single inclusive sample. Representative Feynman diagrams are shown below, in Fig. 

• The loop-induced gluon fusion pp —>■ HZ{ll) -I- 0,1 jets at leading order merged into a single inclusive sample. 
Cf. Fig.j^below for a selection of contributing Feynman diagrams and our definition of the 1-jet contribution 
in this channel. 

Detailed predictions are presented for the invisible Z{ll)H{inv) and hadronic Z{ll)H{bb) Higgs boson decays. 
Using this framework, we also present a realistic phenomenological analysis deriving anticipated LHC Run H 
constraints to the (kj, Ky) coupling parameters. 

This paper is organised as follows. In Section |TTj we outline the basic structures of the Higgs-Strahlung process 
and point out the impact of higher jet multiplicities accounted for through multi-jet merging in a large variety of 
relevant distributions. In Section |HI[ we use our toolkit to explore in detail possible new physics contributions. 
There, we derive the LHC Run H constraints. We draw our conclusions in Section [rv] 

II. ZH PRODUCTION IN THE STANDARD MODEL 
A. Higher-order corrections, multi-jet merging, and simulation set-up 

In this section we discuss in some detail the two dominant ZH production channels, namely quark-induced 
Drell-Yan (DY) type Higgs-Strahlung and the loop-induced gluon fusion (GF), depicted at leading order (LO) in 
Fig. [3 In particular, we study the impact of multi-jet merging to NLO accuracy for the DY and to LO accuracy 
for the GF contributions. 

The DY component comprises the zero- and one-jet squared amplitudes at NLO, illustrated in Fig. [^upper and 
lower panel, respectively. While the zero-jet GF part (Fig. |3(a,b)) is loop-induced and therefore formally LO, in 
the context of multi-jet merging it is more convenient to regard it as a subset of NNLO corrections in the sense 
of counting coupling constant powers wrt. the DY part. In our definition of GF with an extra jet, besides the 
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Figure 2: Upper panel: representative Feynman diagrams for Drell-Yan ZH production at tree level (a) and at one-loop 
level (b). Lower panel: the same for ZHj at tree level (c), one-loop level (d,e), and corresponding real corrections (f). 





Figure 3: Upper panel: representative Feynman diagrams for the Loop^ contribution. While the gluon fusion contributions 
to ZH (a,b) and ZHj (c,d) are indisputably purely loop induced, the squared loop amplitude of diagrams with external 
quarks and a closed fermion loop (e-h) constitute a finite and gauge invariant subset of NNLO corrections to ZHj. The 
latter diagrams of course also interfere with the tree level amplitude and are therefore included, on the amplitude level, in 
the NLO corrections as well (cf. Fig. He)). 


loop-induced diagrams with three external gluons (Fig. Hc,d)), we include all diagrams with a closed quark loop 
and an external quark line (Fig. He-h)). This definition, like the 0-jet gluon fusion component, forms a finite and 
gauge invariant subset of NNLO corrections to ZHj and captures all diagrams which contain a squared Yukawa 
coupling at the squared amplitude level at NNLO QCD. Note that at the amplitude level there is an overlap 
of Feynman diagrams between DY and OF. E.g., diagram Fig. [^e), interfered with the tree level amplitude, 
contributes to NLO DY ZHj, while the same diagram Fig. ^g) is also part of the GF amplitude, contributing at 
loop-squared NNLO. 

Assuming that the invisible sector couples to the Higgs boson only, there is no interference between signal 
and background amplitudes in ZH, H —>■ inv. This is not true for H ^ bb decays, in which case additional 
contributions must be considered. Besides the Higgs decay, bb pairs can be produced through QCD and through 
weak interactions, for example via Z —>■ bb. Accordingly, when the H ^ bb decay is treated as a part of the 
matrix elements as shown in Fig. |^b,c), the amplitude interferes with the tree-level QCD continuum Hl~bb 
production Fig. j^a). Analogously, the tree-loop interference with diagrams of the kind Fig. Qd-g) occurs as a 
background. In order to capture spin correlations and off-shell effects in the gluon fusion ZZ background, we take 
the loop-squared amplitude of diagrams like Fig. Hf,g) into account with the full final state. At this point it is 
worth mentioning that due to spin considerations Z ^ bb and H ^ bb diagrams do not interfere. Some of these 
contributions have not been considered before in the literature. 

While multi-jet merged predictions for the DY channel at NLO have been discussed in [T^] and the merging 
in the loop-induced channel has been technically introduced in m, here we are mostly interested in using this 
technology for detailed studies. In this context, it is worth pointing out that the theoretical precision for the 
DY channel at fixed order is known up to NNLO in the QCD and up to NLO in the electroweak perturbative 
series For the GF contribution, only estimates [16] of NLO corrections in the infinite top mass limit exist; 
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Figure 4: Representative Feynman diagrams which contribnte to the interference of the QCD continnum (a) with the signal 
(b,c) and background (d-h) when the Higgs decays into bb, or when the bb pair is produced via a Z boson. The loop- 
induced ZZgf contribution (f,g) is furthermore included with full l'^l~bb final state in order to capture spin correlations 
and off-shell effects. 


this is due to the fact that a full calculation is hampered by the presence of many scales and, correspondingly, 
a prohibitive complexity in the necessary multi-loop integrals. Similarly to the gluon fusion process for Higgs 
boson and Higgs pair production, the approximation underlying the estimates mentioned above results in a large 
correction factor K ^ 2 to the overall cross section. 

The Sherpa event generator m is used throughout this letter, supplemented with OpenLoops [TB] for the 
calculation of all loop contributions and Collier m for the evaluation of tensor integrals. Finite width effects 
and spin correlations from the leptonic Z boson decay are fully accounted for in the simulation. For the multi-jet 
merging at leading order we employ the ideas of [BD] , adapted for loop-induced processes in . For the merging 
of next-to-leading order matrix elements, the MEPS@NLO algorithm is used. This is implemented along the 
standard LO multi-jet merging algorithms in Sherpa, which also provides tree-level amplitudes, tools for infrared 
subtraction in the calculation of NLO QCD cross sections, and the simulation of parton showers, hadronisation, 
hadron decays, etc. [53]. Throughout our studies we use the NNPDF3.0 parton distribution functions at NLO 
accuracy |24j . In our estimates of theoretical uncertainties, we focus on the usual renormalisation and factorisation 
scale variations by up to a factor of two in the fixed-order part of the simulation. From all these possible scale 
choices we omit the ones in which the factorisation and renormalisation scale prefactors differ by a factor of 4. Since 
there is no higher-order calculation for the GF contribution, we not only consider the usual scale uncertainties, but 
we also consider the effect of higher-orders on the total cross section (ie., the GF rates account for K = 2). Since 
the iF-factor estimate most likely is too naive an assumption, we estimate the associated uncertainty by varying 
the ilT-factor in the range from 1.0 to 4.0. We compare the resulting error bands with the ones obtained from 
the customary scale variations in Section [ll B| In addition, there are two more sources of uncertainties stemming 
from the combination of the fixed-order matrix elements with the parton shower. The first one is related to the 
jet cut used in the multi-jet merging, Qcut, which we vary according to Qcut/GeV € {15, 20, 25}. In addition, the 
uncertainty related to the resummation performed numerically in the parton shower is estimated by varying the 
starting scales with factors in {^0.5,1.0, ^2.0}. 


B. Invisible decays: Z{ll)H{inv) 

Invisible Higgs decays occur in many models collectively referred to as “Higgs-portal” models, see e.g. |5^. In 
these models, the Higgs boson is the mediator between the SM particles and an unknown sector with no other 
tangible interactions with the other Standard Model fields and therefore a prime candidate for Dark Matter. Higgs 
boson production in association with a, Z is a particularly suitable channel for invisible Higgs decay searches due 
to its clean signature with large amounts of missing energy from the undetectable Higgs decay products recoiling 
against the (boosted) leptonic Z boson decays. 

We start the study of the signal sample by applying some typical basic selection cuts. We require two same- 
flavour opposite charged leptons with transverse momentum pr/ > 20 GeV in the pseudo-rapidity range |? 7 i| < 2.5 
and an invariant mass in the region \mii — mz\ < 15 GeV. Jets are reconstructed with the anti-IcT algorithm with 
resolution parameter R = 0.4 and pxj > 30 GeV in the pseudo-rapidity range \yj \ < 5, using FastJet [3S]. After 
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Figure 5: Relative size of the gluon-fusion contribution to the cross section as a function of a minimum cut. We show 

individual curves for the total inclusive cross section and the 1-jet inclusive cross section. Uncertainty bands are obtained 
from scale variations in the gluon fusion sample, keeping the denominator fixed. The NLO Drell-Yan and the loop-induced 
gluon fusion samples are both merged up to one jet, respectively denoted as MEPS@NLO and MEPS@LoOP^ 


applying these kinematic selections, the total signal cross section is strongly dominated by the DY component 
while the GF mode contributes with only 0(10%) to the total ZH cross section [8]. This finding seems to allow 
ignoring the GF channel for all practical purposes. However, selection cuts in searches for invisible Higgs decays 
typically include a minimum requirement, drastically reducing the overwhelmingly large tt-|-jets and U-|-jets 

backgrounds. As shown in Fig. applying such a cut substantially changes the relative composition of the signal 
cross section, and it enhances the relative contribution of the gluon fusion production mode to up to 0(30%) of 
the DY mode. The origin of this increase can be mainly traced back to a harder transverse momentum spectrum 
triggered by a top quark threshold enhancement around mzH ~ 2TOt. This is supported by the finding in the 
left panel of Fig. where we have varied the mass of the heavy quark running in the loop. In contrast, the DY 
contributions do not feature such an enhancement but rather show the typical s-channel suppression for large 
energies. Therefore, despite its small contribution to the inclusive cross-section, the GF mode can become a 
significant player in the boosted regime and a proper modelling of this component is of vital importance. In the 
right panel of Fig. we show the theory uncertainties stemming from scale and AT-factor variations, as detailed 
above, for the DY and the GF mode. Clearly, the AT-factor variation leads to large effects, as the factor of two 
applied in both directions directly translates into an uncertainty, which is about twice as large as the effect of 
the standard scale variation in the GF mode. This size, about 30% or so, is typical for a merged sample at LO, 
especially in view of the fact that it is at least of order C>(a|). In contrast, the scale uncertainty on the DY 
sample is much smaller, about 10-20%. 

In the lower panels of figure (right panel), we also compare the missing transverse energy spectrum of the 
gluon fusion component obtained from a simple LO matrix element plus parton shower simulation (LoOP^-|-PS) 
with the one obtained from a merged calculation (MEPS@LoOP^), taking into account matrix elements with 
up to one extra jet. In the boosted regime above A™'®® > m^, the LoOP^-|-PS simulation significantly under¬ 
shoots the spectrum of the one using merging technology, with the discrepancy reaching around 100% around 
A™'®® ~ 500 GeV and further increasing with energy. This discrepancy has a origin similar to the finite top-mass 
effects in A-|-jets production studied in [10] . namely that the extra jet emissions significantly impact the loop 
structure of the matrix elements. One can artificially suppress this effect by increasing the value of the top quark 
mass, thereby pushing the relevant scale for any loop structure effects to higher energies. This is shown in the left 
panel of figure ^ where the discrepancy between the LoOP^-|-PS and MEPS@Loop^ simulations becomes much 
less severe. 

Apparently, the effects induced by higher multiplicity jet-emission matrix elements are significant and beyond 
the scope of conventional parton showers alone. They can be accounted for by applying matrix element merging 
techniques as demonstrated here, since they correctly fill those phase space regions that are typically problematic 
for the parton shower. This provides a very robust handle on theory uncertainties related to the application of 
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Figure 6: Missing transverse energy distributions after basic selection cuts in the signal channel of the invisible Higgs decay 
search, assnming a branching fraction for H —>■ inv of 1. Left panel: Varying the mass of the top-quark running in the 
loop in the GF contribution shows the threshold effect, extending to the tail of the distribution. Right panel: Uncertainty 
bands obtained from scale variations along with iV-factor variations for the DY and GF contributions. The NLO Drell- 
Yan and the loop-induced gluon fusion are both merged up to one jet. The bottom panel presents the ratio between the 
MEPSOLoop^ to the Loop^-fPS and the MEPSONLO to the MC@NLO. 


vetoes in searches for invisibly decaying Higgs bosons in ZH production, which quite often are an important 
feature in the search strategies. Such vetoes on extra jet emissions are commonly used, cf. for example [B], to 
suppress the backgrounds associated with Higgs-Strahlung, such as top-pair production and similar. Following 
our discussion until now, we anticipate that a jet veto will further suppress the fraction of the loop-induced signal 
component, even when large FI™*®® is being required. A nice way to have some idea about the impact of a jet veto 
is to remind ourselves that the no-emission probability of an additional parton or jet with a transverse momentum 
p_L off a quark g or a gluon g can be roughly estimated, to leading logarithmic precision, using Sudakov form 
factors. Schematically they are given by 


^ 9.3 (mq, P-l) = exp 




-Cf,. 


dgi as{q\) Pq 

^ log “T “ T'9.9 

91 tt \ 


where Cp = 4/3 and Ca = 3 are the colour charges of the quark and gluon and are given by 
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with np the number of active flavours. The occurrence of the colour factors easily motivates why the probability 
for not emitting a jet is larger for quark than for gluon induced processes. 

Defining jet veto efficiencies as 


e(pi) 


^incl ’ 


(4) 


that is the fraction of the inclusive cross section which survives a jet veto applied to jets above a certain transverse 
momentum cut we confront in Fig. [^the simple Sudakov approximation for jet vetoes in the production of 
colour singlet systems, i.e. ZH final states with an invariant mass of m = niH + wz, with the exact results 
stemming from our more detailed simulation. It is remarkable in how far the simple approximation is able to 
reproduce the more exact result in the limit of small transverse momentum cuts applied in the jet veto. The 
results shown in the figure confirm that in the experimentally relevant ranges around 20 GeV, the gluon fusion 
contribution is largely suppressed by jet vetoes due the initial state gluon’s propensity to radiate. 
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Figure 7: Jet veto efficiencies for Drell-Yan-like contributions simulated via MEPS@NLO and for the gluon fusion com¬ 
ponent calculated via MEPS@Loop^. We compare these predictions to simple Sudakov approximations. 


A somewhat alternative way to suppress backgrounds is based on the observation that at large A™'®®, the Z 
and H bosons tend to be more or less back-to-back, rendering a selection cut on their relative azimuthal angle 
ill™®®) an efficient means to improve the signal-to-background ratio. This is because additional jets would 
decorrelate the Z and the H in the transverse plane will effectively be vetoed by such a cut. As shown in Fig.[^ 
these assumptions hold for the Drell-Yan-like contributions. The distributions peak strongly at large azimuthal 
separations A™'®®) and at small values of In the case of the gluon fusion contribution, 

however, the enhancement of the distributions is much less pronounced in these regions, due to the larger level 
of QCD radiation decorrelating the ZH pair; configurations in which the Higgs recoils against a jet rather than 
the Z have a strong impact here. This is even the case at — Pt\/Pt — ^ A(/)(ll, A™®®) = 27r. As can 

be seen in the lower panels of Fig.[^ even in this kinematic regime there is a significant contribution from one-jet 
events in case of the gluon fusion component, whereas this region is depleted of 1-jet events for DY topologies. 
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Figure 8: Differential distributions for the major kinematic background suppression cuts in invisible Higgs decay searches. 
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Performing a search binned in jet multiplicities, as it was done in reference [S], therefore retains sensitivity to the 
gluon fusion component. In fact, as shown in Tab. |lj the gluon fusion component can be as large as 40% after 
applying typical selection cuts in the 1-jet inclusive bin. Modelling this very contribution reliably requires the 
additional jet emission matrix elements and makes the merging techniques applied here an indispensable tool. 


MEPS@NLO MEPSQLoop^ 



^ incl [fb] 






\mii - mz\ < 15GeV,pT! > 20GeV, \yi\ < 2.5 

34.5t?;) 

21.11^3 

13.4l^i 

4.91^4 

l-74t°;«i 

q Q + 1.6 

'^•^-0.9 

> 120 GeV 

9.7ljj 

^•yo_0.69 

4.74l°;^g® 

2 

^•^-0.8 

0 Q5+0.45 

U.i:^0_o.28 

1 Q5+0.97 

-‘-•^'^-0.56 

A0(//, £;??"=) > 2.5 

s.oljj 

4 Q7+0-88 

-0.69 

o q^+0.61 
'5.U^_0.57 


0.95t°:« 

1 42+°-'^4 

-‘-•^^-0.41 

\pt{11) - < 0.25 

6.5lj'2 

/I SI +0-83 

4-o1-0.65 

32 

1 ^7+0.78 
-0.46 

0.881°:" 

0 70+0.37 


Table I: Cut flow for typical selection cuts in invisible Higgs decay searches. We list the individual con¬ 
tributions from Drell-Yan production modes MEPS@NLO and the loop-induced MEPS@LOOP^ component in 
pp —> (H —>■ mv){Z —>■ e^e~, 11 ^pT) at the LHC y/s = 13 TeV. Uncertainties are obtained from scale variations as de¬ 
scribed in the text. For the loop induced contributions, they become as large as the ones one would obtain from varying 
the K-Factor in some cases. This is despite the fact that the if-factor variation error bands of the differential distributions 
in Figs. 1^ and [^exceed the scale variation error bands considerably. 

Apparently however, the gluon fusion component was accounted for neither by ATLAS nor by CMS in their 
searches for invisible Higgs boson decays at Run I 016]. We find that for the ATLAS analysis, the gluon fusion 
component can indeed be neglected. A jet-veto requirement in conjunction with selection cuts similar to the 
ones shown in table |l] suppress the GF component to merely 0(4%) of the total signal rate. The CMS analysis, 
however, takes the 1-jet exclusive bin into account separately, thereby retaining sensitivity to the gluon fusion 
component. We explicitly checked that the significance of the GF contributions in the 1-jet exclusive bin at Run 
I energies is comparable to our findings in table |T| 


C. Hadronic decays: Z{ll)H{bb) 

We now analyse the Higgs-Strahlung channel for the H —>■ bb decay mode. In this case, the Higgs candidate is 
part of a multi-jet system that should contain not only its decay products, but also the associated final state QCD 
radiation. This simple picture is blurred by initial state QCD radiation and additional particles originating from 
the underlying event “splashing” into the fat-jet system stemming from the Higgs boson decay. This complicated 
final state renders proper modelling of the QCD emissions an indispensable requirement for a successful and robust 
analysis of this process. In this section, we will therefore discuss, in particular, the relevance of multi-jet merging 
techniques. 

To highlight the effect of higher order effects and to quantify the impact of multi-jet merging, we follow the 
BDRS analysis [3| as a well understood benchmark: First, we impose some basic selection cuts requiring two 
same-flavour opposite charged leptons with transverse momentum ppi > 30 GeV, pseudo-rapidity \r]i\ < 2.5 and 
invariant mass in the window 75 GeV < mu < 105 GeV. We then impose the reconstructed Z boson to have 
a large transverse momentum pT,a > 200 GeV. In the BDRS algorithm the hadronic final states of events are 
clustered into fat jets using the Cambridge-Aachen algorithm with radius R = 1.2. The analysis demands at least 
one fat jet with ppj > 200 GeV and |r?j| < 2.5, acting as the candidate for the Higgs boson. This candidate is 
tagged through jet substructure techniques including the mass drop criteria and the requirement of three filtered 
subjets for which the two hardest ones need to be 6-tagged. Our analysis assume a flat 70% 6-tagging efficiency 
and a 1% mistag rate. 

In Fig. 0 (left panel) we display the filtered Higgs jet transverse momentum In analogy to the A™'®® in 

the case of invisible decays above, MEPS@Loop^ presents an enhancement with respect to the ZH LoOP^-|-PS. 
This effect is noticeably smaller than in the invisible scenario, however it is still relevant. In the invisible search 
the requirement of A™*®® > 120 GeV leads to smaller invariant masses for the combined ZH system, which in turn 
set the scale for the parton shower to populate the phase space with the emission of extra jets, while here the cut 
on the Z boson transverse momentum is at 200 GeV. The emission phase space offered to the parton shower thus 
is larger in the H > 66 case. As a consequence, in 77 —>■ inv, the exact matrix elements, only included through 
the merging, have a larger impact. Additionally to this kinematic effect, the hadronic Higgs decay is naturally 
more sensitive towards QCD radiation which also induces some differences [36]. 

In Fig. we present some of the QCD radiation dynamics in the invariant mass distribution for the filtered 
Higgs. We first notice that the corrections to the production of the Higgs boson, that are performed here via 
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Figure 9: Transverse momentum distribution of the filtered Higgs candidate (left) and azimuthal angle Hbvks) 

(right) between the Higgs boson and the Z candidate for pp —> {H —^ bb){Z —>■ e^e~ pT) production at y/s = 13 TeV. 
The NLO-merged Drell-Yan contribution is shown in black and the loop-induced LO-merged gluon fusion mode in red. The 
bottom panel presents the ratio between the jet-vetoed and the inclusive sample and the central panel (left only) presents 
the ratio between MEPS@Loop^ merged up to one jet and the parton shower sample. 


the MEPS@NLO and MEPS@LoOP^, contribute to the high mass tail > 125 GeV. Extra parton 

emissions from the fixed order or from the parton shower arising from the ZH production can be reclustered in 
the Higgs fat-jet, therefore enhancing its mass. On the other hand, the invariant mass region < 125 GeV 



Figure 10: Invariant mass distribution of the filtered Higgs fat jet candidate in pp —>■ {H —> bb){Z —>■ e'^e~,p^p~) 

at ^/s = 13 TeV. The NLO DY contribution is shown in black and the loop-induced GF mode in red. Both samples are 
merged up to one jet. The central panel presents the ratio between the MEPS@Loop^ to the MEPS@NLO and the bottom 
panel the ratio between the jet-veto to the inclusive samples. 
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is populated from the shower radiation off the Higgs decays, therefore decreasing the reconstructed invariant mass 
j^BDRS^ Based on these observations we can for instance understand the profile of the gluon fusion contribution 
with respect to the Drell-Yan-like shown in the bottom panel. While their effects at < 125 GeV have the 

same source (i.e. shower emissions off the bb pair), at > 125 GeV we observe a big enhancement that goes 

from 0(18%) to 0(40%) of the DY rate. This is a side effect of the larger radiation pattern arising from the gluon 
fusion component, which benefits in particular from the larger initial state colour charge CaICf = 9/4. Indeed, 
that captures the size of the enhancement in a very good approximation 0.18 x 9/4 ^ 0.4. Besides, we notice 
that the jet veto suppresses the cross-section by approximately a constant factor over the full distribution. 

This clearly indicates that both the parton shower off the Higgs decays and the corrections to the production are 
properly covering all the important phase space regions. 

The i?BDRs) distribution shows that typical selection cuts of order A(j){ll, i?BDRs) > 2.5 have a subleading 

impact on the inclusive rates for the signal component, see Fig. ^ (right panel). Additionally, in the presence 
of an extra jet veto this azimuthal correlation requirement can be pushed even further, almost without extra 
losses. In Table |ll| we display the impact of the cut — mnl < 10 GeV and the extra-jet veto requirement. 

While the cut H^drs) > 2.5 has an imperceptible impact to both components, the extra-jet veto weakens 

the gluon-fusion signal contribution to a subleading level. Thus, if possible one should use other handles (than 
extra-jet vetoes) to the background suppression, especially for possible BSM studies intrinsically associated to the 
loop-induced component. We will further comment on this in the following section. 



MEPSQNLO 

MEPS@Loop2 

cuts 

^incl 

? 

O 

b 

^incl 

1 

o 

b 

BDRS reconstruction 

0.37 

0.18 

0.07 

0.02 

-mnl < 10 GeV 

0.16 

0.09 

0.03 

0.01 


Table II: Cross sections for the Drell-Yan and loop-induced components of pp —> {H —> bb){Z —>■ pApT) production 

at LHC ^/s = 13 TeV. Both samples are merged up to one jet and the selection cuts follow the BDRS analysis that is 
described in the text. The rates are given in fb and account to 70% b-tagging efficiency. Hadronisation and underlying 
event effects are accounted for. 


III. HIGGS-STRAHLUNG: BOOSTING COUPLING CONSTRAINTS 

In this section we analyse the constraining power to the top-quark Yukawa coupling that can be derived from 
thepp —>■ (iJ —>■ hb){Z —>• e+e“,p+p“) production. Following Eq. ^ we notice that the Drell-Yan component does 
not present any dependence on the top Yukawa Kt but only on the size of the HVV coupling kv- On the other 
hand, the gluon fusion develops a dependence on both ny and Kt*. Hence, we can write the matrix element for 
the ZH production as 


M = KtMt + KyMv , 


(5) 


where the t stands for the top Yukawa contributions and V for the contributions proportional to the HVV 
coefficient. The dependence on these coefficients can be straightforwardly translated to the Higgs pth spectrum 
via 


da 

dpTH 


datt 

- - 1 - KtKv 

dpTH 


datv 

dpTH 


+ Kv 


dayy 

dpTH 


( 6 ) 


Therefore, the pth distribution encodes the information about both the size and sign of the top Yukawa. To 
estimate the LHC sensitivity towards these coefficients, we consider the major backgrounds for Higgs-Strahlung, 
namely ti+jets, Zbb+jets, and ZZew- Besides these standard contributions, we also accounted for the loop- 
induced gluon fusion ZZcf production depicted Fig. The interferences with the QCD continuum as described 
in section [IlA[ see Fig. were shown to be subleading in this analysis. 

As for the signal, it is important to properly model the QCD radiation for the background components. The tt is 
generated with the 0-jet bin at NLO and merged up to 3 jets simulated at LO. The Zbb and ZZew are generated 


Only the relative sign between Ky and Kt is physical, thus only positive K,y is considered without loss of generality. 
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pBDRS [QeV] 


Figure 11: Transverse momentum distribution of the filtered Higgs candidate for the Higgs-Strahlung signal and 

backgrounds at ^/s = 13 TeV. We also display the loop-induced gluon fusion component for the BSM hypothesis with 
negative top Yukawa, Kt = —1. 


via MC@NLO and the loop-induced ZZgf at LO. The cut-flow is presented in Tab. |III| We avoid applying extra 
jet veto requirements, since it would deplete the GF signal component as derived in the previous section. Notice 
that the loop-induced ZZqf production presents a non-negligible rate after the BDRS reconstruction, however 
the Higgs mass window selection — mnl < 10 GeV yields a subleading size to it. 

In Fig. E] we present the signal and background transverse momentum distributions after the selection cuts 
depicted in Table III The background components are under control through the whole spectrum with the S/B 
ratio increasing towards higher energies. The negative top Yukawa Kt = —1 displays an amount of events that 



20 40 60 80 100 120 140 

L[fb-'] 


Figure 12: Confidence level for disentangling the negative top Yukawa hypothesis Kt = —1 from the Standard Model. We 
display the results for pp —>■ {H —>■ bb){Z —^ e~, pZp~) based on the distribution. 
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cuts 

ZHgF ACt = — 1 

ZHgf 

ZHdy 

tt-l-jets 

Z6fe-|-jets 

ZZew 

ZZgf 

BDRS reconstruction 

1.48 

0.07 

0.37 

0.29 

13.83 

0.79 

0.10 

- m/f| < 10 GeV 

0.63 

0.03 

0.16 

0.02 

0.35 

0.02 

0.002 


Table III: Cut flow for Z/f+jets (gluon fusion and Drell-Yan-like components), tt+jets, Z66+jets and ZZ+jets. Further¬ 
more, we generate the EW and loop induced QCD components. All simulations were performed with Sherpa-|-OpenLoops. 
The rates are given in fb and account to 70% b-tagging efficiency. Hadronisation and underlying event effects taken into 
account. 


surpasses the Drell-Yan and background components. It largely benefits from the atv term that in the SM 
represents a destructive interference in the whole distribution with a bigger magnitude than the other two 

terms <Jtt and uvv separately. 

In Fig. [^we show the projection of the reach in the ZH analysis for the coupling determination. We analyse 
the information from the different bins via the CLs method [35j and estimate the integrated luminosity 

necessary to exclude the negative top Yukawa solution at 95% CL. A conservative systematic uncertainty of 50% 
is inferred to the GF channel. The BSM hypothesis can be excluded with ~ 30 fb~^ of data. 


IV. SUMMARY 

We have studied the Higgs-Strahlung process at the LHC merging the zero and one jet multiplicities for the 
Drell-Yan and loop-induced gluon fusion via the MEPS@NLO and MEPS@LoOP^ algorithms, respectively. We 
have shown that the multi-jet merging is a fundamental ingredient to properly model the gluon fusion component. 
The merging leads to significant contributions with respect to LO-l-PS simulations. For instance, for typical 
H —>• invisible searches at pth ~ 500 GeV the correction factor is of order 0{2). 

A proper modelling of extra QCD emissions becomes even more important for the H ^ bb decay, since the 
Higgs candidate is part of this multi-jet system. We scrutinised the signal contributions at the boosted kinematics 
and showed that MEPS@NLO and MEPS@LoOP^ provide a good description for the relevant distributions. 
In particular, we observed significant improvements to the transverse momentum p®^^® and reconstructed mass 
for the Higgs candidate. 

Higgs-Strahlung search strategies often rely on extra jet veto requirements that, however, challenge the stability 
of perturbative expansions. We show that MEPS@NLO and MEPS@LoOP^ considerably decrease the impact of 
jet vetoes on the uncertainties in comparison to MC@NLO and LOOP^+PS, respectively. Furthermore, a larger 
suppression to the loop-induced component is observed. At the boosted regime, in particular, extra jet vetoes can 
deplete this signal component to subleading levels. 

Finally, we also estimate the constraining power to the top Yukawa coupling via jet substructure techniques for 
pp —>■ (iL — bb){Z —>■ e+e“, production. We perform a full analysis accounting for all the major background 

components that include, for instance, the EW and loop-induced QCD ZZ production. We conclude that the 
Higgs-Strahlung can be used to access both the size and sign of the top Yukawa coupling. Including conservative 
systematic uncertainties, the Run H LHC can exclude at 95% CL the negative top Yukawa solution Kt = —1 with 
only ~ 30 fb“^. 
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